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ABSTRACT 

We present deep CCD photometry of the globular cluster M13 (NGC6205) 
in the Stromgren uvby(3 system, and determine a foreground reddening of 
E(b — y) = 0.015 ± 0.01 mag. From a fit to the main-sequence of metal-poor 
subdwarfs with Hipparcos parallaxes, we derive (m — M) = 14.38 ± 0.10 which 
implies an age near 12 Gyr assuming [Fe/H] = —1.61 and [ot/Fe] = 0.3. The 
distance-independent ((& — y)o , Co ) diagram indicates that M13 and metal-poor 
field subdwarfs of similar metallicity must be coeval to within ± 1 Gyr. In 
addition, we find that, at any given {b — y) color, there is a large spread in the 
Co index for M13 red-giant branch (RGB) stars. We suspect that this scatter, 
which extends at least as faint as the base of the RGB, is most likely due to 
star-to-star variations in the atmospheric abundances of the CNO elements. We 
also note the existence of what appears to be two separate stellar populations 
on the HB of M13. Among other possibilities, it could arise as the result of 
differences in the extent to which deep mixing occurs in the precursor red giants. 
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1. Introduction 

Despite a concerted effort by many researchers over the past few decades (see, e.g., 
the reviews by VandenBerg, Bolte & Stetson 1996, [Stetson, VandenBerg fc Bolte 19TJS] , 



Sarajedini, Chaboyer fc Demarque 19971 and references therein), it has not been possible to 
reach a consensus on either the absolute or the relative ages of the Galactic globular clusters 
(GCs). Much of the controversy in this field can be traced to the difficulty of determining 
reliable values for the basic parameters that enter into a determination of the cluster age 
such as distance, reddening, overall metallicity and detailed abundance patterns. 

In this respect, the Stromgren uvbyf3 photometric system offers unique advantages 
over broad-band photometry since it can provide precise estimates of T e ff , [Fe/H] and 
surface gravity for F and G stars, as well as the reddening, on a star-by-star basis (see 
Schuster & and Nissen 1989a and references therein). Furthermore, databases of quite 
homogeneous Stromgren photometry exist for large samples of stars, both metal-poor and 
metal-rich (Schuster & Nissen 1988; pisen 1983| , [19841) . Many of these st ars have new, 



high-precision determinations of their parallaxes from the Hipparcos mission, making them 
particularly useful for distance determinations via the main-sequence fitting technique. 
Moreover, since the Stromgren system includes a u filter, which is entirely on the short 
wavelength side of the Balmer jump, it can be used to probe the temperatures and gravities 
of horizontal-branch (HB) stars, including those on extended blue HBs ( |de Boer, Schmidt 
fc Heber 1995ft . 



As is well known, considerable evidence has accumulated during the past ~ 20 years 
which indicates that red giants in globular clusters have mixed much deeper into their 
nuclear-burning interiors than canonical models predict — see, e.g., Kraft et al. (1998), 
Kraft et al. ( 1997}] and the reviews by |Da Costa (1997)| and [Kraft (1994)| for discussions of 



the extensive literature on this subject. Indeed, in order to explain the observed abundance 
patterns of C, N, O, Na, Al, and Mg (also see IBilachowski et al. 199(j|; [Shetrone 199tj|), it 



seems to be necessary for the mixing to penetrate so closely to the H-burning shell that 
some helium must also be dredged up into the surface layers during the evolution (e.g., 
Langer, Hoffman fc Zaidins 19971 , Sweigart 1997a,b, and |Cavallo, Sweigart fc Bell 1998|) . 



As discussed in [Sweigart (1997b)] , even small amounts of He mixing could have large effects 
on the HB and account for at least part of the second-parameter phenomenon (i.e., the 
diversity in HB morphology among GCs having very similar chemical compositions). 
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2. Observations and Data Reduction 

All of the M13 observations were obtained using the Nordic Optical Telescope on La 
Palma, Canary Islands, and Stromgren uvhyf3 filters during the period from June 27 to July 
02, 1995. The weather was excellent during the observing run, giving seeing values between 
0'/43 and 0'/8, with a mean value close to C'60. The detector used was a SITe 1024 square 
chip providing a total field of 3 x 3 arcminutes. All photometric reductions of the cluster 
frames were carried out using DAOPHOT, ALLSTAR, ALLFRAME, and D AO GROW (see 
Stetson 1987, 1990, 1994). In each frame, approximately 20000 stars were detected. For the 
standard stars the residuals in the transformation from the instrumental to the standard 
system were: V: m 004; (b - y): 0^005; m x : m 008; c x : m 010 and (3: 0^008. For M13 the 
error in the zero-point for each filter is estimated to be of order CPOl. The data reduction 
procedures will be described in more detail in a forthcoming paper. 



3. Cluster Parameters 

Before attempting to determine the distance to M13, it is necessary to know the amount 
of interstellar absorption towards the cluster as well as the stellar iron content (preferably 
the detailed chemical abundance patterns). We have used our Stromgren photometry and 
the calibrations of Schuster & Nissen (1989a, Eq. 1) to determine the foreground reddening. 
A number of their calibrating stars were included as standards in our observing program to 
ensure that our photometric system would be close to theirs. The M13 stars used included 
turn-off and subgiants with ^-magnitudes between 17.5 and 20.5 (~ 500 stars). In this 
way, we derive a reddening of E(b — y) = 0^015 ± 0.01, or E(B — V) — 0.021, which is quite 
typical of the values found for M13 (e.g, [Harris 1996| ). 



In principle, the cluster metallicity can also be determined from Stromgren photometry. 
Using (again) the calibrations of |Schustcr fc Nissen (1989a) , we derive a mean metallicity 



of [Fe/H] = — 1.7 for ~ 500 stars. However, we found some indication of a possible error 
in the transformation of our instrumental v magnitudes Q to the standard system as a 
function of mi and, since an error of 0^01 in the mi index corresponds to an error of ~ 0.1 
dex in the estimated metallicity, the [Fe/H] value that we have derived for M13 cannot be 
considered robust. As a result, we have adopted [Fe/H] = — 1.6 ± 0.2 as a compromise of 
our photometric determination and recent spectroscopic results ( |Kraft et al. 1998] , Carretta 



fe Cratton 1997] ) . Finally, we adopt [a/Fe] = 0.3 which seems to be appropriate for the halo 



5 No such trends were evident for the other filters or indices. 
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population ( |Carney 1996| ). There are no measurements of this quantity in main-sequence 
or turnoff stars in Ml 3, but see [Kraft et al. [T997] and [Pilachowski fc Armandroff (1996) 
for pertinent discussions. 



3.1. The Cluster Distance and Age 

Ideally, the distance to M13 would be determined by a main-sequence fit of the C-M 
diagram for the unevolved cluster stars to a large sample of local subdwarfs having the 
same metallicity as the cluster and very precise parallaxes. Unfortunately, the number of 
subdwarfs with [Fe/H] ~ —1.6 and with relative errors in their parallaxes at the few percent 
level remains small. Potentially-useful subdwarf candidates were selected from the list of 
[Nissen, H_0g fc Schuster (1997)| , with metallicities derived from Stromgren photometry and 
parallaxes determined by Hipparcos . As in the case of M13, the individual reddenings and 
metallicities of the subdwarfs were determined using the calibrations of (Schuster fc Nissen 



1989a 



Following the usual procedure (see, e.g., |VandenBerg et al. 1996| ), the dereddened colors 
of the subdwarfs were individually adjusted to compensate for the (small) differences in 
their metallicities relative to that of M13 using the isochrones of [VandenBerg et al. (1998 



as transformed to the Stromgren system using [Kurucz ( 1992 )| color-temperature relations. 



In the fit of the M13 ((it — y)o, V) and ((u — y)o, u) C-M diagrams to the resultant "mono- 
metallicity" subdwarf sequence, we used stars with cr(7r)/7r < 0.07, where 7r is the parallax, 
and photometric metallicities in the range —1.69 < [Fe/H] < —1.57. Consequently, the 
model-based color corrections are all quite small for these stars. Furthermore, by adopting 
such stringent constraints on the relative parallax error, the Lutz-Kelker corrections ( [Lutz 



|fc Kelker 1973| , |Lutz 1979] ) amount to only a few hundredths of a magnitude, at most, and 



possible systematic biases for the whole sample are negligible as well (see the discussion by 
Pont et al. 1998) |. We applied Lutz-Kelker corrections to the subdwarfs using the formula 



(> 



Two of the stars used in the main-sequence fit (HD 64090 and HD 188510) are listed by Gratton et 



al. (1997) and Pont et al. (1998~)| as spectroscopic binaries or radial-velocity variables. However, Carney 



ct al. (1994) has 19 velocity measurements for HD 188510 over 8.2 years, and they find no indication of 
radial-velocity variations (rms. 0.8 km/s), while five out of six measurements of the velocity of HD 64090 
by |5tryker et al. (1985)| are in excellent agreement. |Stryker et al. (1985)| quotes a semiamplitude limit of 



14 km/s for their sensitivity, their one discrepant velocity value differs from the mean of the other five by 
27.4 km/s. Moreover, regardless of which C-M diagram is used in the subdwarf fit, these two stars show no 
sign of being "overluminous" with respect to the other two that we have used (HD 34328 and HD 25329) and 
the MS fiducial. Therefore, if HD 64090 and/or HD 188510 truly are binaries, the secondary component(s) 
must be much fainter than the primaries. 
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[] given by |Hanson (1979) with n = 2. Figure 1 illustrates the ((it — y)o, V) diagram for 



M13 with the subdwarfs overplotted. The main-sequence fit yields (m — M)y = 14.44 and 
(m — M) u = 14.49 with an uncertainty of ±0.10 mag. 

Isochrones from [VandenBerg et al. ( 1998)1 f° r P e /H] = —1.61, [a/Fe] =0.3, and ages 
of 10, 12, and 14 Gyr have also been superposed on the data in Fig. 1. These indicate that 
the age of M13 is close to 12 Gyr. The reduction in age from values near 14-16 Gyr that 
were commonly obtained a few years ago is due less to revisions in the cluster distance than 
to improvements in the stellar models. Indeed, the M13 distance modulus derived here is 
within 0.1 mag of the values determined by e.g., |Richer fc Fahlman (1986)| and [Buckley "fc" 



Longmore (1992)] in the pre-Hipparcos era. It is, instead, the use of improved opacities 



and an equation of state that treats non-ideal effects (e.g., |Profht 1993[ ), the adoption of a 



revised bolometric correction scale (see |VandenBerg 1997| ), and the assumption of [a/Fe] 



= 0.3, instead of a scaled-solar heavy-element mix, that is primarily responsible for the 
reduced GC ages. It should be noted that the present models do not include the effects of 
helium diffusion, which can be expected to reduce the age at a given turnoff luminosity by 
~ 1 Gyr (see, e.g., Proffit fc VandenBerg 1991). 



4. The ((& — y) , c ) Diagram 

Having extensive Stromgren photometry for M13 stars fainter than the main-sequence 
turnoff allows us to determine whether the cluster differs significantly in age from field stars 
of similar metallicity. For F and G stars, Cq =C\ — 0.2E(b — y) is a measure of the size of 
the Balmer jump, which, in turn, is a measure of surface gravity and hence of evolutionary 
state. Therefore, the c\ index can be used to distinguish between evolved and unevolved 
stars and to determine the ages of stars near the turnoff (see, e.g., Fig. 7 in |VandenBcrg 
fc Bell 1985| ) independently of any knowledge of their distances. This approach has been 
used by |Schuster fc Nissen (1989b")| to determine ages for a large sample of metal-poor 
field stars. In Figure 2, we have overplotted the photometry for stars in their catalogs 
having metallicities between —1.4 and —1.9 on the ((b — y)o , Co ) diagram for M13. It is 
immediately obvious from this figure that the field stars follow the cluster locus closely, 
indicating that the two stellar populations must have very nearly the same age. The 
maximum difference between the turnoff Cq values is ±0™025, which translates into an age 
difference of ~ ±1 Gyr, which, we emphasize, is completely independent of the adopted 



7 Because of the high precision of the parallaxes for the selected stars, the choice of n has no significant 
impact on the derived distance modulus. 
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distance scale. We cannot, of course, exclude the possibility that some of the field stars on 
the lower main sequence are older, only that none of the evolved stars has a greater age 
than M13. 

A surprising feature of the data plotted in Fig. 2 is that, at any given (6 — y)o color, 
the M13 giant-branch stars exhibit a large spread (CP1 - CP15) in the Cq index. As we 
have no spectra for the observed stars, it is not clear at this point how these Cq variations 
should be interpreted. However, as noted in §1, it is well known that stars on the M13 RGB 
encompass large variations in the C, N, O, Mg, and Al abundances, most likely due to deep 
mixing during the post-main-sequence evolution. Star-to-star variations in these elements 
are the most likely cause of the Co variations, in view of the fact that several investigations 
( [Bond 1980| , [Anthony-'l'warog fc Twarog 1994| , [Anthony-'l'warog, Twarog fc Craig 1995| ) 



have demonstrated that stars with strong CH and CN bands can have significantly lower 
Ci indices than similar stars of normal CH and CN strength. Spectroscopic observations of 
the stars showing extreme c\ values are needed to identify the elements responsible for the 
variations. 

If the large spread in the c\ observations is due to variations in the elemental 
abundances due to mixing, then Fig. 2 would suggest that such mixing commences near the 
base of the red-giant branch. One cannot help but speculate that a thermal instability in 
the H-burning shell might be at least partly responsible for the observed scatter since |Von 



[Rudloff, VandenBerg fc Hartwick (1988)| have shown that the greatest potential for such an 
instability occurs just as a star begins to ascend the giant branch. This possibility should 
be investigated further. 



5. The Horizontal Branch 

We now turn to a brief discussion of the properties of the observed horizontal branch. 
In Fig. 1, a ZAHB locus (VandenBerg et al. 1998), which represents the extension of 
the isochrones to the core He-burning phase, has been superposed on the data on the 
assumption of the adopted distance modulus and reddening. It is apparent that the coolest 
of the HB stars are reasonably well-matched by the model ZAHB. However, at (u — y) 
< 0.95, the HB stars appear to be ~ (KM "overluminous" compared to the theoretical 
locus. We note that this "anomaly" has also been seen by Dorman et al. (1998, private 
communication) in their HST observations of M13 (using the F336W and F555W filters). 
In fact, a preliminary analysis of Stromgren data that we have in hand for NGC288 and 
NGC 6752, which both possess blue HBs, show the same unexpected morphology. This 
cannot be explained in terms of evolution away from a ZAHB since the evolution is slowest 
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near the zero-age locus and most stars should, therefore, be found adjacent to it. Unless 
the discrepancy is due to a problem with the color-T e g- relations derived from model 
atmospheres, which seems unlikely, a different ZAHB must apply to these stars; i.e., M13 
(as well as NGC 288 and NGC 6752) must possess at least two distinct HB populations. 

At first sight, these observations appear to be in remarkable agreement with the 
recent predictions by Sweigart (1997a,b), who has demonstrated that, if red-giants undergo 
sufficiently deep mixing as to modify their envelope helium abundances appreciably, then 
their decendants on the HB should be bluer and more luminous than predicted by canonical 
models. In fact, there is some indication in his calculations (compare, e.g., Figs. 7 and 9 in 
Sweigart 1997b) that, for the hottest HB stars, the He-mixed models and those based on 
standard assumptions, will nearly coincide. Thus, the near match of the canonical ZAHB 
with the lower-bound of the M13 stars fainter than V ~ 17.7 (see Fig. 1) can be understood. 



On the other hand, |Rood et al. (1998)1 have found that there is no apparent discrepancy 
between canonical HB models and M13's HB when comparisons between the two are made 
on C-M diagrams involving bands that are bluer than the HST F336W filter (which is close 
to the Stromgren it). So, at the present moment, it is not clear how this will be resolved. 

Finally, we note that in a (u — y, u) diagram both the ZAHB appropriate to the coolest 
BHB stars and the locus of subgiant stars are very nearly horizontal and very well defined; 
thus it is easy to measure the u luminosity difference {Au) to high accuracy. From the 
isochrones of [VandenBerg et al. (1998)1 we estimate that Au changes at a rate of ~ 0.065 



mag/ Gyr. Given the high precision available here we estimate that this method will provide 
the possibility to measure the age difference between clusters of similar metallicity to better 
than 0.5 Gyr. 



6. Summary 

The capability of Stromgren photometry to reveal important new features of GC C-M 
diagrams is undeniable. In this investigation of M13, we have found that: 1) The reddening 
of M13 is E(b - y) = 0.015 ± 0.01, which corresponds to E(B - V) = 0.021 ± 0.014. 2) 
The cluster distance is (m — M) = 14.38 ± 0.10, based on C-M diagram fits to local 
subdwarfs with Hipparcos parallaxes. This implies an age of 12 ± 1.5 Gyr according to 
recent theoretical isochrones for [Fe/H] = —1.61 and [a/Fe] = 0.3 (but which do not take 
He diffusion into account). 3) There is a large spread in c% for stars from the base to the tip 
of the RGB, which is probably indicative of star-to-star variations in the CNO abundances. 
4) Field stars with metallicities similar to that of M13 are found to be coeval with M13 to 
within pa ±1 Gyr. 5) The HB shows evidence for two distinct populations, which may be 
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related to the observed scatter in ci on the RGB. It is also possible that we are seeing the 
signature of He mixing in upper RGB stars, as discussed by Sweigart (1997a,b). 
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Fig. 1.— The (u — y)o, V) color-magnitude diagram of M13, on which 10, 12, 14 Gyr 
isochrones (for [Fe/H] = —1.61 and [a/Fe] = 0.3), as well as local subdwarfs having Hipparcos 
parallaxes (filled circles), have been superposed. In the order of increasing My, the subdwarfs 
are: HD 34328, HD 64090, HD 188510, and HD 25329. (Adding in HD 19945 and HD 29907 
(open squares), which have [Fe/H] = —1.91 and [Fe/H] = —1.63, respectively results in the 
same distance modulus). The isochrone/ZAHB colors were adjusted by 5(u — y) = 0.05 mag 
in order to achieve coincidence of the predicted and observed lower main sequences. We 
attribute the failure of the models to match the RGB to problems with the Kurucz (1992) 
color transformations: R. Bell (1998, private communication) has obtained much better 
agreement of the same isochrones with the M13 C-M diagram (in the vicinity of the turnoff) 
using his color-T eff relations. 

Fig. 2. — The ((b — y) , cq ) diagram of M13, with field stars having —1.9 < [Fe/H] < —1.4 
from [Schuster fc JNissen (19895"J| overplotted as filled circles. Note the large scatter in cq on 



the cluster RGB (but not among the few field giants that are present in the plot) at any 
given (b - y) . 
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